(* Content-type: application/vnd.wolfram.mathematica *)

(*** Wolfram Notebook File ***)
(* http://www.wolfram.com/nb *)

(* CreatedBy='Mathematica 8.0' *)

(*CacheID: 234*)
(* Internal cache information:
NotebookFileLineBreakTest
NotebookFileLineBreakTest
NotebookDataPosition[       157,          7]
NotebookDataLength[     20666,        502]
NotebookOptionsPosition[     19874,        472]
NotebookOutlinePosition[     20212,        487]
CellTagsIndexPosition[     20169,        484]
WindowFrame->Normal*)

(* Beginning of Notebook Content *)
Notebook[{

Cell[CellGroupData[{
Cell[BoxData[{
 RowBox[{
  RowBox[{"\[CapitalOmega]", "=", 
   RowBox[{"4", "*", 
    RowBox[{"\[Pi]", "/", "3"}]}]}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"MAXN", " ", "=", " ", "15"}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"pos", " ", "=", " ", 
   RowBox[{"Array", "[", 
    RowBox[{
     RowBox[{"\[Theta]", "&"}], ",", "MAXN"}], "]"}]}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{"pos", "[", 
    RowBox[{"[", "1", "]"}], "]"}], " ", "=", " ", "0"}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{"For", "[", 
  RowBox[{
   RowBox[{"i", "=", "1"}], ",", 
   RowBox[{"i", "<", " ", "MAXN"}], ",", 
   RowBox[{"i", "++"}], ",", "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"d", " ", "=", " ", 
     RowBox[{"Array", "[", 
      RowBox[{
       RowBox[{
        RowBox[{"1", "/", "2"}], "&"}], ",", 
       RowBox[{"i", "+", "1"}]}], "]"}]}], ";", "\[IndentingNewLine]", 
    RowBox[{
     RowBox[{"d", "[", 
      RowBox[{"[", 
       RowBox[{"i", "+", "1"}], "]"}], "]"}], "=", "\[Delta]"}], ";", 
    "\[IndentingNewLine]", "\[IndentingNewLine]", 
    RowBox[{"u", " ", "=", " ", 
     RowBox[{"Array", "[", 
      RowBox[{
       RowBox[{
        RowBox[{
         RowBox[{"-", "1"}], "/", 
         RowBox[{"Log", "[", 
          RowBox[{"\[Epsilon]", " ", 
           RowBox[{"d", "[", 
            RowBox[{"[", "#", "]"}], "]"}]}], "]"}]}], "&"}], ",", 
       RowBox[{"i", "+", "1"}]}], "]"}]}], ";", "\[IndentingNewLine]", 
    RowBox[{"u1", " ", "=", " ", 
     RowBox[{"Array", "[", 
      RowBox[{
       RowBox[{
        RowBox[{"u", "[", 
         RowBox[{"[", "#", "]"}], "]"}], "&"}], ",", "i"}], "]"}]}], ";", 
    "\[IndentingNewLine]", "\[IndentingNewLine]", 
    RowBox[{"G", " ", "=", " ", 
     RowBox[{"DiagonalMatrix", "[", 
      RowBox[{"Array", "[", 
       RowBox[{
        RowBox[{
         RowBox[{"1", "/", 
          RowBox[{"(", 
           RowBox[{"8", " ", "\[Pi]"}], ")"}]}], "&"}], ",", 
        RowBox[{"i", "+", "1"}]}], "]"}], "]"}]}], ";", "\[IndentingNewLine]", 
    RowBox[{"For", "[", 
     RowBox[{
      RowBox[{"j", "=", "1"}], ",", 
      RowBox[{"j", "\[LessEqual]", 
       RowBox[{"i", "+", "1"}]}], ",", 
      RowBox[{"j", "++"}], ",", "\[IndentingNewLine]", 
      RowBox[{"For", "[", 
       RowBox[{
        RowBox[{"k", "=", 
         RowBox[{"j", "+", "1"}]}], ",", 
        RowBox[{"k", "\[LessEqual]", 
         RowBox[{"i", "+", "1"}]}], ",", 
        RowBox[{"k", "++"}], ",", "\[IndentingNewLine]", 
        RowBox[{
         RowBox[{
          RowBox[{
           RowBox[{"G", "[", 
            RowBox[{"[", "j", "]"}], "]"}], "[", 
           RowBox[{"[", "k", "]"}], "]"}], " ", "=", " ", 
          RowBox[{
           RowBox[{
            RowBox[{
             RowBox[{"-", "1"}], "/", "\[Pi]"}], "*", 
            RowBox[{"Log", "[", 
             RowBox[{"Sqrt", "[", 
              RowBox[{"2", "-", 
               RowBox[{"2", " ", 
                RowBox[{"Cos", "[", 
                 RowBox[{
                  RowBox[{"pos", "[", 
                   RowBox[{"[", "k", "]"}], "]"}], "-", 
                  RowBox[{"pos", "[", 
                   RowBox[{"[", "j", "]"}], "]"}]}], "]"}]}]}], "]"}], 
             "]"}]}], " ", "+", " ", 
           RowBox[{"1", "/", 
            RowBox[{"(", 
             RowBox[{"8", "*", "Pi"}], ")"}]}]}]}], ";", 
         "\[IndentingNewLine]", 
         RowBox[{
          RowBox[{
           RowBox[{"G", "[", 
            RowBox[{"[", "k", "]"}], "]"}], "[", 
           RowBox[{"[", "j", "]"}], "]"}], " ", "=", " ", 
          RowBox[{
           RowBox[{"G", "[", 
            RowBox[{"[", "j", "]"}], "]"}], "[", 
           RowBox[{"[", "k", "]"}], "]"}]}], ";"}]}], "\[IndentingNewLine]", 
       "]"}]}], "\[IndentingNewLine]", "]"}], ";", "\[IndentingNewLine]", 
    "\[IndentingNewLine]", "\[IndentingNewLine]", 
    RowBox[{"G1", " ", "=", " ", 
     RowBox[{"Array", "[", 
      RowBox[{
       RowBox[{
        RowBox[{
         RowBox[{"G", "[", 
          RowBox[{"[", "#1", "]"}], "]"}], "[", 
         RowBox[{"[", "#2", "]"}], "]"}], "&"}], ",", 
       RowBox[{"{", 
        RowBox[{"i", ",", "i"}], "}"}]}], "]"}]}], ";", "\[IndentingNewLine]",
     "\[IndentingNewLine]", 
    RowBox[{"p", " ", "=", " ", 
     RowBox[{"u", ".", "G", ".", "u"}]}], ";", "\[IndentingNewLine]", 
    RowBox[{"p1", " ", "=", " ", 
     RowBox[{"u1", ".", "G1", ".", "u1"}]}], ";", "\[IndentingNewLine]", 
    "\[IndentingNewLine]", 
    RowBox[{"ubar", " ", "=", " ", 
     RowBox[{
      RowBox[{"Total", "[", "u", "]"}], "/", 
      RowBox[{"(", 
       RowBox[{"i", "+", "1"}], ")"}]}]}], ";", "\[IndentingNewLine]", 
    RowBox[{"u1bar", " ", "=", " ", 
     RowBox[{
      RowBox[{"Total", "[", "u1", "]"}], "/", "i"}]}], ";", 
    "\[IndentingNewLine]", "\[IndentingNewLine]", 
    RowBox[{"v", " ", "=", " ", 
     RowBox[{
      RowBox[{"\[CapitalOmega]", "/", 
       RowBox[{"(", 
        RowBox[{
         RowBox[{"(", 
          RowBox[{"i", "+", "1"}], ")"}], "*", "\[Pi]", "*", "ubar"}], 
        ")"}]}], " ", "+", " ", 
      RowBox[{
       RowBox[{"\[CapitalOmega]", "/", 
        RowBox[{"(", 
         RowBox[{
          RowBox[{
           RowBox[{"(", 
            RowBox[{"i", "+", "1"}], ")"}], "^", "2"}], "*", 
          RowBox[{"ubar", "^", "2"}]}], ")"}]}], "*", "p"}]}]}], ";", 
    "\[IndentingNewLine]", 
    RowBox[{"v1", " ", "=", " ", 
     RowBox[{
      RowBox[{"\[CapitalOmega]", "/", 
       RowBox[{"(", 
        RowBox[{"i", "*", "\[Pi]", "*", "u1bar"}], ")"}]}], " ", "+", " ", 
      RowBox[{
       RowBox[{"\[CapitalOmega]", "/", 
        RowBox[{"(", 
         RowBox[{
          RowBox[{"i", "^", "2"}], "*", 
          RowBox[{"u1bar", "^", "2"}]}], ")"}]}], "*", "p1"}]}]}], ";", 
    "\[IndentingNewLine]", "\[IndentingNewLine]", 
    RowBox[{"F", " ", "=", " ", 
     RowBox[{"Limit", "[", 
      RowBox[{
       RowBox[{
        RowBox[{"(", 
         RowBox[{"v", "-", "v1"}], ")"}], "*", 
        RowBox[{"Log", "[", "\[Delta]", "]"}]}], ",", 
       RowBox[{"\[Delta]", "\[Rule]", "0"}], ",", 
       RowBox[{"Direction", "\[Rule]", 
        RowBox[{"-", "1"}]}], ",", 
       RowBox[{"Assumptions", "\[Rule]", 
        RowBox[{"\[Epsilon]", ">", "0"}]}]}], "]"}]}], ";", 
    "\[IndentingNewLine]", 
    RowBox[{"k", " ", "=", " ", 
     RowBox[{"Maximize", "[", 
      RowBox[{
       RowBox[{"{", 
        RowBox[{"FullSimplify", "[", 
         RowBox[{"Re", "[", 
          RowBox[{"F", "/.", 
           RowBox[{"\[Epsilon]", "\[Rule]", "0.1"}]}], "]"}], "]"}], "}"}], 
       ",", "\[Theta]"}], "]"}]}], ";", "\[IndentingNewLine]", 
    "\[IndentingNewLine]", 
    RowBox[{
     RowBox[{"pos", "[", 
      RowBox[{"[", 
       RowBox[{"i", "+", "1"}], "]"}], "]"}], " ", "=", " ", 
     RowBox[{"\[Theta]", "/.", 
      RowBox[{"k", "[", 
       RowBox[{"[", "2", "]"}], "]"}]}]}], ";"}]}], "\[IndentingNewLine]", 
  "]"}]}], "Input",
 CellChangeTimes->{{3.518440783317552*^9, 3.518440783877502*^9}, {
   3.518440814706459*^9, 3.518440882250953*^9}, {3.518440987500345*^9, 
   3.518441035654086*^9}, {3.518441113753427*^9, 3.518441115975266*^9}, {
   3.518441149394973*^9, 3.518441221926009*^9}, {3.518441255470563*^9, 
   3.518441373680977*^9}, {3.518441410968585*^9, 3.518441648293933*^9}, {
   3.518445825039301*^9, 3.51844591366432*^9}, {3.518449787493126*^9, 
   3.518449807593542*^9}, {3.518449838582437*^9, 3.518449845009269*^9}, {
   3.518449911437523*^9, 3.518449912637981*^9}, {3.518450012165515*^9, 
   3.518450027357726*^9}, {3.518450213196362*^9, 3.518450245189562*^9}, {
   3.518450304072066*^9, 3.51845030552689*^9}, {3.518450409252792*^9, 
   3.518450521305021*^9}, {3.518450940012388*^9, 3.518450946762911*^9}, {
   3.518451018335815*^9, 3.51845107056742*^9}, {3.518451130298596*^9, 
   3.518451147479322*^9}, {3.518451229442954*^9, 3.518451232855136*^9}, {
   3.518451402935868*^9, 3.518451425363926*^9}, {3.518451558379286*^9, 
   3.518451560630882*^9}, {3.518451654633902*^9, 3.518451672111032*^9}, {
   3.518451727144075*^9, 3.518451727398678*^9}, {3.518451985887266*^9, 
   3.518451993815848*^9}, {3.518452111876958*^9, 3.518452159462903*^9}, {
   3.518452524980255*^9, 3.518452526874404*^9}, 3.518452787630545*^9, {
   3.518453015946513*^9, 3.518453016907118*^9}, {3.518453548992123*^9, 
   3.518453551506025*^9}}],

Cell[BoxData[
 RowBox[{
  StyleBox[
   RowBox[{"Maximize", "::", "objv"}], "MessageName"], 
  RowBox[{
  ":", " "}], "\<\"The objective function \[NoBreak]\\!\\(\[Infinity]\\)\
\[NoBreak] contains a nonconstant expression \[NoBreak]\\!\\(\[Infinity]\\)\
\[NoBreak] independent of variables \[NoBreak]\\!\\({\[Theta]}\\)\[NoBreak]. \
\\!\\(\\*ButtonBox[\\\"\[RightSkeleton]\\\", ButtonStyle->\\\"Link\\\", \
ButtonFrame->None, ButtonData:>\\\"paclet:ref/Maximize\\\", ButtonNote -> \
\\\"Maximize::objv\\\"]\\)\"\>"}]], "Message", "MSG",
 CellChangeTimes->{3.518452672621449*^9, 3.518452922095433*^9, 
  3.518453151553087*^9, 3.51845370487534*^9}],

Cell[BoxData[
 RowBox[{
  StyleBox[
   RowBox[{"ReplaceAll", "::", "reps"}], "MessageName"], 
  RowBox[{
  ":", " "}], "\<\"\[NoBreak]\\!\\({\[Theta]}\\)\[NoBreak] is neither a list \
of replacement rules nor a valid dispatch table, and so cannot be used for \
replacing. \\!\\(\\*ButtonBox[\\\"\[RightSkeleton]\\\", ButtonStyle->\\\"Link\
\\\", ButtonFrame->None, ButtonData:>\\\"paclet:ref/message/ReplaceAll/reps\\\
\", ButtonNote -> \\\"ReplaceAll::reps\\\"]\\)\"\>"}]], "Message", "MSG",
 CellChangeTimes->{3.518452672621449*^9, 3.518452922095433*^9, 
  3.518453151553087*^9, 3.518453704879287*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[{"pos", "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"posplot", " ", "=", " ", 
   RowBox[{"Array", "[", 
    RowBox[{
     RowBox[{
      RowBox[{"{", 
       RowBox[{
        RowBox[{"pos", "[", 
         RowBox[{"[", "#", "]"}], "]"}], ",", "1"}], "}"}], "&"}], ",", 
     RowBox[{"Length", "[", "pos", "]"}]}], "]"}]}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{"ListPolarPlot", "[", "posplot", 
  "]"}], "\[IndentingNewLine]"}], "Input",
 CellChangeTimes->{{3.518441037789459*^9, 3.518441085483683*^9}, {
   3.518441122970692*^9, 3.518441124405947*^9}, {3.518441184123223*^9, 
   3.518441225369895*^9}, {3.518441264124657*^9, 3.518441286465093*^9}, 
   3.518441318430005*^9, {3.518441494487455*^9, 3.518441494875827*^9}, 
   3.518449849878002*^9, {3.518449921453828*^9, 3.518449957286841*^9}, {
   3.5184500313082*^9, 3.518450097983718*^9}, {3.518450139637193*^9, 
   3.51845018893172*^9}, {3.518450273248931*^9, 3.518450300702606*^9}, {
   3.518450527073717*^9, 3.518450569560631*^9}, {3.518450637228312*^9, 
   3.51845064813587*^9}, {3.518450690385146*^9, 3.518450690860974*^9}, {
   3.518450807963967*^9, 3.518450842145078*^9}, {3.518450930502316*^9, 
   3.518451002640815*^9}, {3.518451042357578*^9, 3.518451063882701*^9}, {
   3.518451211343136*^9, 3.518451214689672*^9}, {3.518451349465169*^9, 
   3.518451350595764*^9}, {3.518451535866664*^9, 3.518451536594248*^9}, {
   3.518452052245864*^9, 3.518452094691427*^9}, {3.518452128131464*^9, 
   3.518452128345092*^9}, 3.518452344475679*^9, {3.518452419904983*^9, 
   3.518452485204054*^9}, 3.518452997188685*^9, {3.51845383802328*^9, 
   3.518453839878891*^9}}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{"0", ",", "3.141592637037476`", ",", "1.570796303257783`", ",", 
   RowBox[{"-", "1.570796352441104`"}], ",", 
   RowBox[{"-", "0.7853981845401939`"}], ",", "0.8955225467338861`", ",", 
   RowBox[{"-", "2.4550031105093173`"}], ",", "0.46130405885271386`", ",", 
   RowBox[{"-", "3.863486870345311`"}], ",", 
   RowBox[{"-", "2.0113748638243605`"}], ",", 
   RowBox[{"-", "0.42397820772344075`"}], ",", "2.0237534886693216`", ",", 
   RowBox[{"-", "3.4665489308620483`"}], ",", 
   RowBox[{"-", "1.1966784297142077`"}], ",", 
   RowBox[{"\[Theta]", "/.", "\[VeryThinSpace]", "\[Theta]"}]}], 
  "}"}]], "Output",
 CellChangeTimes->{
  3.518441225691422*^9, {3.518441264851208*^9, 3.518441286754397*^9}, {
   3.518441318830057*^9, 3.518441331919357*^9}, {3.518441496039103*^9, 
   3.51844151089238*^9}, {3.518445856115266*^9, 3.518445879966301*^9}, 
   3.51844985022451*^9, {3.518449928452064*^9, 3.518449957863287*^9}, {
   3.518450031786667*^9, 3.518450085927179*^9}, {3.51845014266699*^9, 
   3.518450189462215*^9}, 3.518450230850542*^9, {3.518450275525293*^9, 
   3.518450301005984*^9}, 3.518450615925774*^9, {3.518450686233732*^9, 
   3.518450691688697*^9}, {3.518450815671942*^9, 3.518450842627622*^9}, {
   3.518450951067722*^9, 3.518451003867343*^9}, {3.518451043961738*^9, 
   3.518451073515602*^9}, 3.518451214988404*^9, 3.518451351279728*^9, 
   3.518451537067796*^9, 3.518451643994793*^9, 3.518451710353913*^9, 
   3.518451824255261*^9, {3.518452045969875*^9, 3.518452095351318*^9}, {
   3.518452128628201*^9, 3.518452153659154*^9}, 3.518452227036933*^9, 
   3.5184523428828*^9, {3.518452444122922*^9, 3.51845248567026*^9}, 
   3.51845270566332*^9, 3.5184529990979*^9, 3.518453840381414*^9}],

Cell[BoxData[
 GraphicsBox[{{{}, {}}, {}, GraphicsComplexBox[CompressedData["
1:eJxTTMoPSmViYGCQAWIQDQEf7GGs/2Dwfj/DAZNTd48G2v2Pl1K86hpp9w8i
bv9//tO0yZ3R+/5C1anUV8TNW/DMvvhHax6Q3m+2Zf7tFwxP7Fc3PfNS+/XC
3mHOm4TT+1/sd7WcFfLe48n+Nf2tk7ctf2Nf6RJ99kTNHfvcvinrelhe7A+b
lPI6VPWp/Q77Q3u0vG/v/yo+d8WWj2/2f/ghYb9W6639qxXu0aYht/bvnPNK
RYjpzv7gnAebnhx9Y7/EsPuJS/C7/Ux9HSEP8q7YG4V0cPYnX7d3kH+4K+bk
2/3D3X8AIWL8lA==
    "], {{}, 
     {Hue[0.67, 0.6, 0.6], 
      PointBox[{15, 16, 17, 18, 19, 20, 21, 22, 23, 24, 25, 26, 27, 
       28}]}, {}}], {}, {{{}, {}, {}, {}}, {}}},
  Axes->True,
  AxesOrigin->{0, 0},
  CoordinatesToolOptions:>{"DisplayFunction" -> ({
      Sqrt[Part[#, 1]^2 + Part[#, 2]^2], 
      Mod[
       ArcTan[
        Part[#, 1], 
        Part[#, 2]], 2 Pi]}& ), "CopiedValueFunction" -> ({
      Sqrt[Part[#, 1]^2 + Part[#, 2]^2], 
      Mod[
       ArcTan[
        Part[#, 1], 
        Part[#, 2]], 2 Pi]}& )},
  DisplayFunction:>Identity,
  PlotRange->{{Automatic, Automatic}, {Automatic, Automatic}},
  PlotRangeClipping->True,
  PlotRangePadding->Scaled[0.02]]], "Output",
 CellChangeTimes->{
  3.518441225691422*^9, {3.518441264851208*^9, 3.518441286754397*^9}, {
   3.518441318830057*^9, 3.518441331919357*^9}, {3.518441496039103*^9, 
   3.51844151089238*^9}, {3.518445856115266*^9, 3.518445879966301*^9}, 
   3.51844985022451*^9, {3.518449928452064*^9, 3.518449957863287*^9}, {
   3.518450031786667*^9, 3.518450085927179*^9}, {3.51845014266699*^9, 
   3.518450189462215*^9}, 3.518450230850542*^9, {3.518450275525293*^9, 
   3.518450301005984*^9}, 3.518450615925774*^9, {3.518450686233732*^9, 
   3.518450691688697*^9}, {3.518450815671942*^9, 3.518450842627622*^9}, {
   3.518450951067722*^9, 3.518451003867343*^9}, {3.518451043961738*^9, 
   3.518451073515602*^9}, 3.518451214988404*^9, 3.518451351279728*^9, 
   3.518451537067796*^9, 3.518451643994793*^9, 3.518451710353913*^9, 
   3.518451824255261*^9, {3.518452045969875*^9, 3.518452095351318*^9}, {
   3.518452128628201*^9, 3.518452153659154*^9}, 3.518452227036933*^9, 
   3.5184523428828*^9, {3.518452444122922*^9, 3.51845248567026*^9}, 
   3.51845270566332*^9, 3.5184529990979*^9, 3.51845384038487*^9},
 ImageCache->GraphicsData["CompressedBitmap", "\<\
eJztW1toFFcY/pKNa6mXQCmWig8B0aalNy1qJSk+ZWN6eWhfarEaNkMgtqYb
kxQqkoJBKU3bVFNTWAUt1YdeCKYvoiJWC9FVi3jplWpbk4og2YbmJkncnc6c
nRNnz545Z5zZxBmch7ns/P/3zX/+/z9n/jNn55Voc11tfbR5gxItqWyMNtRt
UJpKIrFG7VKoACio07ZoCfRzFaC7Q1BVFSGQH9qZvh0kezgT7SeHYqipIfx6
9AtsfuYhLHi9HV09/ZjQJGNXd6D8iTU4OZCGersfx1vXo6GtEx/GqrH9ZBK3
NZ0ia/oD5DBX2w/h7JuLsWJvH6ElahMD+P1SH0bTKobP1+PJ5Z34a1zF+J87
sOypTbgwogqYdxuiW1c+w6ZaBYqioJDsQ1DqtqD7xjjRKCT7ItaCELk8U9uP
4cpHz2LRunMY0u8wmMAbC5fg06tj8nYJWEdxqaEUpbGLGNGRIxfw9mOPo/ny
qGFQCGmt2WfP9BH272yyjuGPNs3W9dm2tmu2UtYziV6seiFOWIv0XRdLPYjE
mkVYnqHONDFs+P9z/E3830H9T513+tS1SVbN5m5z7swgufPzwQ68F12L+sM3
TTkR1nImiRPbomhuj+OTphp88MO/upwS60441XONEGvn3W5S2Z+ilUw6/Hju
HzTJ1Ism46xvGqxRjqAxtI+gwUmc7tXPDcRGjsGV1m3eyDc4IkcwBhMEATB6
jJkC5nsd60BkJTLH1Ut23SeiJuOauU+vlKnn9Onn5QimT5sRBED1mD4tYPaI
B6dKtIQTmVaZek5ktsoRTGRsIJgYCRDBwOAH0eRUwKLMvzMVSGH4cieUl2vQ
0hHHx+/UYNuRm6SuF5TMlN6i1DdNB8Zx/ZtGtHyfRErVdXbiucXVSAwK6eOG
KHc+MCtrPpDZrIv+MO1w/YdQPn8p9vROiG77JTmIK/5M/R5G6r+fcPjigObM
FJLHqrGwtAHnhzM1+Aydu4vltK73KedE39dY27oL8V3b8W6sFd/+NoK0Kpo4
7TeQFiV/yECqyRNQ2n7BLcMUCZtFnU/NJKaXPo0VZWUoKytHVawHA2mRmR7p
FB4WddJMNT0aHpGp5zwa5skRzKPBjCAAqsc8EATMHvFgILISBY9dn4uCAPpc
NKUBfJAcil299Y1QDhfl4mqOiS/mgbciD+3LvKLUqkf94PgNNvWSm6pX5CU3
vBF3TTMcZJSitO4O5fcdfLnR9qko63l+rbJhE48r4ryJWYnmaEmhwoaOZBIi
9oVjWqlbZOskcrdYz5EqbOhI5lGyFHFIayNbxAs9JrfczXpORCzmzf1eshx8
7oqmypG9hebHlaOFpohYLJmfWo+8ThkrnDelcIqKAWeiWXQQNM0+v+KoH+Do
PZChIOOo/gRBjb4ryM5sZs2xmMM9N8/X5kzDPabjWr7bMZsfkpnZAWPWcoOA
3fuAMSExAsa722uca6+60AtP+608Mi76XxS8OvC5KAigz0U2A/iwITKXVvtk
6jkLAHvlCGYBwAaCWQoQILzkd4+I/N9/OziZOV+mnpOZj8oRTGbaQDCZKUC4
8oq/F9PNK9DBQnmwUH6/i/w/IrsSLaNd3TScvy9TzxnOW+QIZji3gWCGcwEi
yAy/imo56VcpU7f4x7kQwfvHuRjB/e95kBm+EL3FSasqmXpOWq1mY031mGTi
LSUwCCaZBAiPeNAjIt9+w2bje6/g27Tg2zSeCAX/A08kiHM=\
\>"]],

Cell[BoxData[""], "Print",
 GeneratedCell->False,
 CellAutoOverwrite->False,
 CellChangeTimes->{{3.518441058114231*^9, 3.518441086088538*^9}, {
  3.518441118447107*^9, 3.518441151600881*^9}, {3.518441184883525*^9, 
  3.518441199381977*^9}}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[""], "Input",
 CellChangeTimes->{{3.518451834562649*^9, 3.518451860813417*^9}, {
   3.518451945288825*^9, 3.518451970527875*^9}, 3.518452002358443*^9}],

Cell[BoxData[
 RowBox[{"Real", "[", 
  RowBox[{"9.357388972221505`", "\[VeryThinSpace]", "+", 
   RowBox[{"1.7661057888493454`", " ", 
    RowBox[{"Log", "[", 
     RowBox[{
      RowBox[{"(", 
       RowBox[{"1.`", "\[VeryThinSpace]", "+", 
        RowBox[{"0.`", " ", "\[ImaginaryI]"}]}], ")"}], "-", 
      RowBox[{"1.`", " ", 
       RowBox[{"Cos", "[", 
        RowBox[{
         RowBox[{"(", 
          RowBox[{"1.`", "\[VeryThinSpace]", "+", 
           RowBox[{"0.`", " ", "\[ImaginaryI]"}]}], ")"}], " ", "\[Theta]"}], 
        "]"}]}]}], "]"}]}], "+", 
   RowBox[{"1.7661057888493454`", " ", 
    RowBox[{"Log", "[", 
     RowBox[{"1.`", "\[VeryThinSpace]", "+", 
      RowBox[{
       RowBox[{"(", 
        RowBox[{"0.9999999999999998`", "\[VeryThinSpace]", "+", 
         RowBox[{"0.`", " ", "\[ImaginaryI]"}]}], ")"}], " ", 
       RowBox[{"Cos", "[", 
        RowBox[{
         RowBox[{"(", 
          RowBox[{"1.`", "\[VeryThinSpace]", "+", 
           RowBox[{"0.`", " ", "\[ImaginaryI]"}]}], ")"}], " ", "\[Theta]"}], 
        "]"}]}], "-", 
      RowBox[{
       RowBox[{"(", 
        RowBox[{"1.8260767243101635`*^-8", "+", 
         RowBox[{"0.`", " ", "\[ImaginaryI]"}]}], ")"}], " ", 
       RowBox[{"Sin", "[", 
        RowBox[{
         RowBox[{"(", 
          RowBox[{"1.`", "\[VeryThinSpace]", "+", 
           RowBox[{"0.`", " ", "\[ImaginaryI]"}]}], ")"}], " ", "\[Theta]"}], 
        "]"}]}]}], "]"}]}]}], "]"}]], "Output",
 CellChangeTimes->{{3.518451835960806*^9, 3.518451861472046*^9}, {
  3.518451948164454*^9, 3.518451970910493*^9}}]
}, Open  ]],

Cell[BoxData[""], "Input",
 CellChangeTimes->{{3.518441126987332*^9, 3.518441128559519*^9}, {
  3.518441191207983*^9, 3.518441191851861*^9}}],

Cell[BoxData[""], "Input",
 CellChangeTimes->{{3.51844120149353*^9, 3.518441201793439*^9}}]
},
WindowSize->{740, 607},
WindowMargins->{{142, Automatic}, {Automatic, 23}},
FrontEndVersion->"8.0 for Linux x86 (32-bit) (February 23, 2011)",
StyleDefinitions->"Default.nb"
]
(* End of Notebook Content *)

(* Internal cache information *)
(*CellTagsOutline
CellTagsIndex->{}
*)
(*CellTagsIndex
CellTagsIndex->{}
*)
(*NotebookFileOutline
Notebook[{
Cell[CellGroupData[{
Cell[579, 22, 8368, 215, 715, "Input"],
Cell[8950, 239, 646, 12, 43, "Message"],
Cell[9599, 253, 603, 11, 43, "Message"]
}, Open  ]],
Cell[CellGroupData[{
Cell[10239, 269, 1634, 30, 88, "Input"],
Cell[11876, 301, 1739, 28, 50, "Output"],
Cell[13618, 331, 3988, 75, 368, 1084, 28, "CachedBoxData", "BoxData", \
"Output"],
Cell[17609, 408, 240, 5, 23, "Print"]
}, Open  ]],
Cell[CellGroupData[{
Cell[17886, 418, 164, 2, 30, "Input"],
Cell[18053, 422, 1567, 40, 55, "Output"]
}, Open  ]],
Cell[19635, 465, 141, 2, 30, "Input"],
Cell[19779, 469, 91, 1, 30, InheritFromParent]
}
]
*)

(* End of internal cache information *)

